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ABSTRACT 

We consider 2-dimensional QCD on a cylinder, where space is a circle. We 
find the ground state of the system in case of massless quarks in a 1/N ex- 
pansion. We find that coupling to fermions nontrivially modifies the large 
N saddle point of the gauge theory due to the phenomenon of 'decompact- 
ification' of eigenvalues of the gauge field. We calculate the vacuum energy 
and the vacuum expectation value of the Wilson loop operator both of which 
show a nontrivial dependence on the number of quarks flavours at the leading 
order in 1/N. 
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Introduction and Summary 



The large N expansion provides a valuable tool for obtaining qualitative in- 
sight into gauge theories [[j]]. The 1/N expansion in gauge theories is based on 
the discovery by 'tHooft that for SU (N) gauge theories with fundamental 
fermions the N dependence of a typical vacuum diagram Eh,l for large iV is 

Eh,l oc N*{±) h {±) l (0.1) 

where H is the number of handles and L the number of holes (fermion loops) 
in the diagram. It is easy to see that the dominant contributions to the 
vacuum energy occurs when H = L = and that it is proportional to iV 2 . 
This simply reflects the fact that there are N 2 gluons in the theory. The 
fermionic contribution to the vacuum energy, on the other hand, must have at 
least one fermion loop, and the leading contribution is _Eo,i oc N. This, again, 
reflects the fact that there are N fermions in a fundamental representation. 
The iV 2 -dependence of the vacuum energy has been explicitly calculated in 
soluble large N matrix models || which, like gauge theories, also possess 
N 2 degrees of freedom. On the other hand, two-dimensional QCD (QCD2) 
on the plane, which is another soluble model at large N |3[], possesses no 
dynamical gauge degrees of freedom and has vacuum energy oc N coming 
from the N fermions. See also the more recent works on pure QCD2 and 
on QCD 2 with fermions || 0. 

In the present paper we continue our study of large N QCD 2 on a 
cylinder. Here space is a circle and hence the gauge field does not decouple. 
After fixing gauge appropriately the gauge field can be described by the N 
eigenvalues of Ai which satisfy fermi statistics. The standard wisdom, based 



on Eqn. ( |0.1| ) and the discussion in the last paragraph, would suggest the 
following two-step procedure for solving the theory at large N — (1) to solve 
the pure gluon theory first (the large iV saddle point of this is easy to de- 
termine and is described by a constant density of eigenvalues; see remarks 
before ( |3.47| ) below) and (2) treat the fermion dynamics subsequently as fluc- 



tuations in the fixed external gauge field background determined in step (1). 
In this paper we explicitly show that such a procedure is incorrect. The rea- 
son for this, basically, is that in presence of quarks there are gauge-invariant 
operators which (unlike the Wilson loop operator) are not periodic functions 
of the eigenvalues and this effectively leads to a noncompact range of the 
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eigenvalues unlike in the pure Yang-Mills theory. This is the phenomenon of 
'decompactification' described in great detail in [0]. It is clear that a constant 
density of eigenvalues is not normalizable in a noncompact domain. Indeed 
we find that the quarks lead to a harmonic oscillator potential for the eigen- 
values, resulting in expressions for physical quantities that are nontrivially 
different from those of pure Yang-Mills theory at the leading N order. We 
present our main results below. 

The vacuum energy in our theory is given by (see Eqn. ( |3.43j ) below) 

E = N 2 [^nj£ + o(l/N)} (0.2) 

where rif is the number of quark flavours and g = g\Z~N as usual denotes 
the scaled coupling constant. The vacuum expectation value of the Wilson 
loop operator W m = (l/N)TrU m , U = exp[i J L A x dx] is given by (see Eqn. 

(W m ) = 2{-l) m (l+—^)J {x m )+o{l/N), x m = 2vr 3 / 4 a4k « = \JgL/27c 

(0.3) 

In the above E = E \/N~ and L = LyN denote an additional iV-scaling 
necessary in order to have a well-defined large N limit in our theory. The 
expression for the eigenvalue density is given by Eqn. (|3.44j ). 

Note the nontrivial dependence of the above expressions on the number 
of quark flavours which clearly shows that adding quarks changes the leading 
large N result. For comparison, note that the vacuum expectation value of 
the Wilson loop operator in pure Yang-Mills theory in two dimensions is (see 
(B below ) (Wm)y M = S mfi + o{l/N). 

This paper is organized as follows. In Sec. 1 we write down the largrangian, 
the path integral and the hamiltonian for QCD on a two-dimensional cylin- 
der. We follow the notation and results of [0. In Sec. 2 we find the ground 
state of the theory in the 1/N expansion and reduce the calculation of the 
vacuum energy to a problem of iV interacting fermions in a harmonic oscil- 
lator external potential. The solution of the latter problem is presented in 
Sec. 3 in the 1/N expansion which allows us to calculate the vacuum energy 
of the full system. We find in the process that we need to scale both the 
radius of the cylinder and the time in a certain way (see Eqns. (|3.29|) and 
(|3.31|) ) to have a well-defined large iV limit. Besides the vacuum energy we 
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also calculate the vacuum expectation value of the Wilson loop operator. In 
Sec. 4 we conclude with some comments about possible relevance to four 
dimensions. 

1 The Action and the Hamiltonian 

We consider the gauge group U(N). As usual we denote the gauge fields, 
which are hermitian matrices, by A" where \x = 0, 1 is the Lorentz index 
and a, b = 1,2, ... ,N are colour indices. The fermions are denoted by ipf a 
where i — 1, 2, . . . , n/ is the flavour index and a = ±1 is the dirac index. We 
consider one space and one time dimension where the space dimension is a 
circle of length L. 

The theory is described by the following path-integral 

Z = j VA (x, t)VA 1 {x, t)Th{>(x, t)Vip\x, t) exp[iS{A , A x , ft)] 

S = tfdtf L dx( - \trF^F^ + ^(iD^ - m^O (1.1) 

F 01 = E = d t A x - d x A + ig[A , A x ], D^ = 3^ + igA^ 

In this paper we will only consider the case of massless quarks, m = 0. 
Let us fix the gauge 

A 1 (x,t) = Dmg[\ a (t)] (1.2) 
The gauge-fixed path-integral becomes [[5] 

Z = J V\(t)Ap(\(0))Ap(\(T))Vilj(x,t)V^(x,t)exp[i(S Y M + S F )} (1.3) 

where 

Sf = Jo dt J Q L dxtyUM + iph 5 (id x - g\)ip] + S coul (1.4) 
Seoul = - Ea.fe Jo dt J L dx J L dyp ab (x)p ba (y)K ab (x - y) 
The kernel K ab (x) is given by (for x G [— L, L]) 

K ab {x) = (^ 2 /4)e^(^-^)[L/(2sin 2 <Xa ~^ ) - txcot n{Xa x ~ Xb) - \x\], a ^ b 
K aa {x) = {g 2 /±)[L/Q-\x\+x 2 /L\ 

(1.5) 
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The action in ( |Q ) is invariant under the "large" gauge transformations M, 



ip a (x,t) -> Q ab (m,x)ip b (x,t), tt ab (m,x) = exp[-igxm a X a ] 

(!-6) 

A (t) A?(t) = A (t) + m a A 

where 

A„ ^ (17) 

The parameter A defines the range [0, A ] of integration of the eigenvalues 
A a in ( |1.3|) . A simple way of understanding the period (|1.7| ) is to consider 
the Wilson loop operators 

W m = ^TrU rn , U = expfa f L Mdx) (1.8) 
iV jo 

which in the gauge ( |1.2| ) evaluate to 

W ™ = T7 X] exp(z^LA a ) ee — ^ exp(i2vrA a /A ) (1.9) 

-'"a a 

Mesonic operators or Hp in (|1 . 1 3 ) below, which involve gluons as well as 
quarks, are not periodic under A a = — > \ a = \o per se but one needs to 
simultaneously transform the quarks also according to ( |1.6|) . This ultimately 
leads to a decompactification of eigenvalues in the effective theory of the 
gluons. 

The factors Ap(A(0)) and Ap(A(T)) in the measure are defined by 

A P (A) ee J] sin[7r(A a - A 6 )/A ] (1.10) 

a<b 

These factors in the measure imply that the initial and final wavefunctions 
are completely antisymmetrized with respect to the A a 's (since to start with 
they must be symmetric with repsect to permutation of the A a 's on account 
of Weyl-symmetry) . This gives rise to the well-known fermionic nature of 
these eigenvalues. 

The hamiltonian corresponding to the action (1^4) is given by 

H = H YM + H F , (1.11) 
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where 



and 



Hp — Hpfl + H cou i 

Hf, = E a Jo dx ^ a j 5 (-id x + g\ a )r (1-13) 
H cou i = - J2ab Jo dx J L dy p ab (x)p ba (y)K ab (x - y) 



The kernels K ab [x — y) are given by ([TT 

The physical Hilbert space of the above theory satisfies the zero charge 
condition (for each colour) 

Q a = J dxp aa (x) = (1.14) 



2 Ground State 



In this section we construct the ground state of the hamiltonian ( |1.11| ) in the 
1/N expansion. We do it in two steps, (a) We first consider the gauge field 
as external and discuss the dynamics of the fermions for fixed A Q 's. We show 
that for any fixed A a 's the dirac sea built out of free fermions is the ground 
state of this problem modulo 1/N 2 corrections, (b) Next we use this result 
to construct the ground state of the full problem where both fermions and 
the eigenvalues are dynamical. 

We should remark that the expectation value of H in the full ground state 
of the theory was already presented in 0. In the discussion below we will 
allow some small overlap with H for the sake of completeness. 



2.1 Fermion dynamics in an external background of 
gauge fields 

We start by discussing eigenstates of Hp t0 (eqn. ( |1.13Q ) in the presence of 
fixed background values of {A a }. We will include the effect of H cou i later on. 
For simplicity we will also work with rif — 1 at first. 
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Since Hp^ is quadratic, the ground state is simply given by filling the 
fermi sea according to the single-particle spectrum 

E a na = 2 -fsgn(a)(n+^) (2.1) 

Eqn. ( p.lD is obtained by noting that the dirac equation [id t + rf{id x — 
g\ a )]ipa = is solved by 

27T7i 

<&t) = £exp[-i££ t + i-— x]C,„ (2-2) 

n L 

with E^ n as in ( [2.1|) . Here we have chosen the convention 7 = <Ji and 
7 1 = —2(72, so that 7 5 = 7°7 1 = 03. a = ±1 represent the right-, left-moving 
fermions respectively. 

In filling the fermi sea, a priori the fermi levels could be different for dif- 
ferent chirality and different colour. However, as we noted in |7J], translation 
invariance of the vacuum demands that the right fermi-momentum must be 
left fermi-momentum minus one for each colour. Let us denote the right 
fermi-momentum for colour a as p F = 2-Kn a F jL. The fermi (dirac) sea is 
therefore the state \np) defined by 



^R,n\n F ) = 


0, n > n a F 




0, n <n a F 


V L , n \n F ) = 


0, n < n a F 




0,n > n a F 


■ defined by 


(|2~^) and sa' 



(2.3) 



relation 

The reader might wonder how the fermi levels can be different for differ- 
ent colours and still produce a colour- singlet state. Indeed the question is 
further complicated by the fact that the fermi levels n F , having come from 
momentum labels of fermions, are non- gauge-invariant; under the large gauge 
transformation ( |1.6| ) they transform as 

n F _> = n a F -m a (2.5) 
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We will defer the detailed discussion of gauge-invariance of the vacuum 
till Sec. 2.2. For the moment, let us check that the state \n F ) satisfies the con- 
straint (P..14Q . Note that Q a = J2a=R,L Qa,a Where Q a ^ a = L Em=-oc ^rn^a-n 

Evaluating this on the state \np) we get 



Qr= E !> Ql= E 1 (2-6) 

m=— oo m=nj,+l 

The sums are divergent. By using the gauge- invariant exponential regulator 
[0] exp[-e|m + A a /A |] we get 

Q a R (e) = - + l+n a F + ^+o(e), Q a R (e) = - - \ - n F - ^ + o(e) (2.7) 
This leads to 

Q a (e) = Q£(e)+Q£(e) = ^ (2.8) 

This is a divergent c-number term which can be removed by normal ordering. 
Thus : Q a := Q a ( e ) — (2/e) = 0. Note that we have been able to achieve this 
without assuming any special values of n F . 

Construction of ground state of H F 

It is easy to see that \n F ) is an eigenstate of the dirac part H F0 of the 
hamiltonian (|1.13|) 

H Ffi \n F ) = E Ffi \n F ) (2.9) 

oo 

£f,o = sAoE E (m + A a /Ao)tr,( 7 5 Q^) (2.10) 

a m=— oo 

where 

Q a m ^ = L(n F \^ m r p Jn F ) = 9(n F - rn) l -±f- + 9(m - n% - 

(2.11) 

Regarding the action of H cou i on \n* F ), it is more convenient to use the 
alternative form 

Hcoui = E H<*0 (2-12) 

a,/3 
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where 

tt L ■s—^ pa,ab,nPf3,ba,-n /„ , „s 

a " ~ 2A§ A [n+(A fl -A b )/A ]* (2J3) 

In the above 

oo 

Pa,ab,n= £ <mC,m+n ( 2 - 14 ) 
m=— oo 

It is a lengthy but straightforward calculation to show that 

[H RR + H LL }\n F ) = E coul \n F ) (2.15) 
where ^ 

C ° Ul ~ 4-Ao £ ^nfe-oo + (A- " A 6 )/Ao] 2 

where the QJi have been introduced above in ( |2.11p . 

The left-right mixing terms H LR = H RL take \n F ) to an orthogonal state 
(which is a linear combination of states with two holes and two particles). If 
we treat H FR as a perturbation term, the first order perturbation correction 
is zero since 

(n F \H LR \n F ) = (2.17) 
Second order perturbation theory gives the contribution 

g ™ m 



^ 2 Ag jx [m + n F -n b F + (A a - A 6 )/A 
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For finite N there is no reason to regard this as a perturbation. However 
in the large N limit we scale the coupling constant g as g = g/VN so that 
A = AoviV. This makes E LR down by 1/N compared to E coul because of 
the two extra powers of Ao in the denominator. This is actually the story 
with conventional scaling of coupling constant. As we shall see in Sec. 3, the 
scaling in our theory also involves scaling of the length L of the circle. This 
in fact brings down Ep R by an additional factor of 1/N. 

Combining the above results, we find that the dirac sea is an eigenstate 
of H F to the leading order in 1/N, with the energy given by 

H F \n F ) = E F \n F ) + o{l/N), E F = E Ffl + E coul (2.19) 
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The sums in (|2.1Q ) and fl2.16| ) are divergent. Using the exponential regulator 



once again as in ( [2.71) we get 



E F , (e) = 9XoJ2i-- e - T§ + + o(e)] (2.20) 



a 



and 



E coul (e) = 4^rE l^f { c n + 2 ( lne " !) + * (^) + 
47TA afc esm ^(<U-<UJ A 

(2.21) 

where 



Keg(^) = (f + (2.22) 

-Kreg( ^\ ^ ) = [2C + ^(^afe) + lp(-W ab ) + W ah {ljj' {w ah ) ~ 1p'(-W ab )}} 
Ao 

(2.23) 

Here u> a & = (£„— £&) / Ao and C is Euler's constant. Also, ij){x) = (d/ dx) In T(x) 
and = 9 x ip(x), T(x) being the standard gamma-function. Note that 

Ep only depends on the gauge-invariant combination 

^ a = \ a + n F X (2.24) 

The appearance of this variable is responsible for decompactification of the 
eigenvalue A a . We will discuss this in more detail shortly. 

A remark is in order here justifying our definition of the regularized quan- 
tities K e g and K reg . The divergent piece in E Ft0 (e) is a constant and, because 
of the constraint that the total number of eigenvalues is N, does not affect 
the dynamics. In the grand canonical ensemble the equivalent statement is 
that such a divergence can be cancelled by a simple additive renormalization 
of the chemical potential term. Similar remarks can be made about the ln(e) 
piece in E cou[ (e) and the constant finite pieces in E F;0 (e) and E coul (e). The 
justification for ignoring the l/(esin 2 ) piece in E cou i(e) is more subtle. The 
main point is that this 1/e comes multiplied with the quadratic pole (£ a — £f>) -2 
which is the only singularity of (|2.21| ) in the limit £ a — > £&. If we take the 
same limit in the unregulated expression ( p. 16 ) we find a quadratic pole with 



residue En(Qn) 2 = E n Qn = Qa- Here we have used (Q a n ) 2 = Q£Va,n. Thus 
the 1/e in (|2.21| ) can be identified with the total charge which must vanish 
by ( |TT1 . 
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Combining all this we get the following regularized expression for the 
eigenvalue Ep 

E F = gX EaiU^o + 1/2) 2 + jfc E a , b K rcg (^) (2.25) 

where K Teg is defined in (|2.23|) . 

In the above we have proved that the filled fermi sea is an eigenstate of Hp 
at leading 1/N order. How does one argue that it is actually the ground state 
of the system? If we ignore for the moment H cou i it is obvious that the filled 
fermi sea \np) is the ground state of H F0 (for any given set of n a F and A a , or 
equivalently, for any given set of £ a 's, which is probably a more appropriate 
specification of external background in our problem). In presence of H cou i the 
argument is not so simple. Let us present several independent reasons why 
the filled fermi sea should be the ground state, (a) In [7) we have shown that 
the expectation value of the meson bilocal operator M xy (X,t) (see eq. (55) 
of [[7J) in this state provides the unique lowest energy translation invariant 
solution to the classical equation of motion, (b) The four-fermi interaction 
represented by H cou i is repulsive in nature, (c) It is easy to show that small 
number of gauge-invariant mesonic excitations always increase the energy of 
this state, (d) In the scaling that we describe in the next section (involving 
g and L) H cou i is subleading to Hp t o by a factor of 1/N. Thus, to leading 
order in 1/N the filled fermi sea must be the ground state. 

2.2 Effective hamiltonian for gauge fields and ground 
state of the full theory 

Now that we have computed the ground state of Hp for fixed external gauge 
field, let us construct the full ground state by the method of separation of 
variables. What we will do now is similar in spirit to solving the Schrodinger 
problem for a central force in three dimensions where we look for solutions 
which are products of radial and angular wave-functions. We solve the an- 
gular problem first and find the eigenfunctions of the angular momentum 
operator. The centrifugal term L 2 /r is evaluated by using the eigenvalue of 
the angular momentum at a fixed r. This is then put back in the full lapla- 
cian to derive an effective hamiltonian for the radial problem which is then 
solved by appropriate methods. 
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The schematic correspondence between the above example and our prob- 
lem will be ip <-> {6,4>), A <-> r, Hf[iP, ^ , A] <-> L 2 (6>, <%, <9</,)/r, \tif) <-> 
^0,0(^5 0) an d Hym ^ —(l/r 2 )d r (r 2 d r ). Note the important fact that |n~p) 
is independent of the A a 's as is clear from the defintion Q2.3|). 

Let us write a product wavefunction of the system as 

|#) = |n» <g> (2.26) 
Using ( ILTTD , flTigp and flT25p we find that 



= |nlr) ® fleff$(A) (2.27) 

where 

^ = ErK)+ £ F ( 2 - 28 ) 

with .Ep given by ( |2.25| ). 

Our first guess at the full ground state would be = \np) ® $o(A) 
where $o(A) is the ground state of if e fj. However, as we have discussed 
in great detail in @ such a state \^) is not gauge- invariant because of the 
shift ( |2.5| ) of the fermi levels under large gauge transformations. The correct 
gauge-invariant ground state is given by a sum over fermi levels 

l*o>=E|n»®*g(X) (2-29) 

Hp 

where the wavefunctions $^(A) are given by 

4°j(A)=« (0) (A + ^) (2.30) 

Here £ G (—00, 00) is the ground state wave-function for the hamil- 

tonian fl2.28| ) where in the kinetic term the operator d\ a has been replaced 
by d^ a . In other words, satisfies the differential equation 

Hu^(0=E u^(0 (2.31) 

where 

W = E^(-3L) + M3 (2-32) 
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In the above Eq denotes the lowest eigenvalue of the hamiltonian. We shall 
discuss in the next section how to evaluate it. 

Note that in expectation values computed in the full theory, the implica- 
tion of the sum-of-product structure (|2.29|) of the full wave-function is that 
the effective range of the eigenvalues becomes the real line. Let us calculate, 
for instance, the expectation value of the full hamiltonian H in the state 
l*o>- 

<*oi# i* > = e n / Ao ^ *S*$tf**S$ = n r « (0) *(o^ (o) (o 

(2.33) 

where H eS and TC are given by ( |2.2£j| ) and ( |2.32| ). The last line shows that 
the presence of fermions forces a decompactification of the eigenvalues A a 
to the gauge-invariant combination A a + n^Ao = £ a which is a real number 
G (—00, 00). It also tells us that there is no fixed fermi level in a compact 
gauge theory with dynamical gauge fields; rather, the effective theory of the 
gauge fields A in the fermi vacuum is given by a density matrix p(np, X\n' F , A'). 
The sum over rip in the above equation corresponds to taking trace over this 
density matrix. 



3 The Large- N Expansion 

In this section we will discuss the large N limit in detail and present the 1/iV 
expansion for some physical quantities. For concreteness, we will consider 
the partition function 

exp(-/3F) = Z = Trexp(-/3if) (3.1) 



where H = Hym + Hp is given by (1.1 1|) . We will also evaluate towards the 



end of this section the vacuum expectation value of the Wilson loop operator 
TiU m defined by flOl). 

The large iV limit involves defining a scaled coupling constant 



9 = gVN (3.2) 

which is held fixed as iV —>■ 00. According to ( |L7| ), we must also define a 
scaled eigenvalue-period 

Ao = Vn\ , Ao = — (3.3) 
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This necessitates scaling of A a , £ a and p a to 

A a = X a /VN, C = U^N, Pa = Pa/VN 



(3.4) 



Note that in terms of the scaled variables p a = —{if 'N)d\ , implying that h 
is 1/JV. Let us, for example, rewrite Eqns. (|2.31| ) and (|2.28|) for the ground 
state energy and eigenf unctions in terms of the scaled variables 



^ (0) ({U)=^ (0) ({C}) 

Ti- = H.ym + 'Hi + %2 
H Y M={f\o/4:Nir)j: a (-d 2 -) 

^i = ^AoEa(UAo + l/2) 2 



(3.5) 



(3.6) 



where the K reg has been defined in ( j2.23| ) (note that in ( [2.23[ ) w a b is also 
equal to [£, -£ b ]/A ). 



3.1 (3 — > oo limit and the ground state energy 

In the (3 — > oo limit, the only state contributing to the trace in (|3.1|) is 
the ground state, discussed in the previous section. Thus the free energy F 
simply coincides with the ground state energy E of (|3.5|) . 

Let us now discuss how to determine Eq. If the two-body interaction Ti.2 
was absent, Eq would be simply given by 

£ ,o = -^iV 3 / 2 (3.7) 

Proof: The ground state for the iV-fermion hamiltonian TCym + ^i is given 
by the Slater determinant 

4(0 = ^D<i^(e i+ i) (s-s) 

where </> n (£),n = 0, . . . , oo are normalized single-particle wavefunctions sat- 
isfying 

d 2 - 

- ^ + 2 *(y Q + \) 2 }M0 = enUO (3-9) 
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It is easy to find explicit solutions of the above equation: 



2 n n\ — 

0„(£) = (^=)- 1/2 H n (V^x) exp(-ux 2 /2), cu 2 = 4tt (3.10) 
ay JSuj 

where the variable x is defined by 

(e/A ) + (1/2) = ax, a = X^N^ 4 (3.11) 

In the above H n are the standard Hermite polynomials. The energy eigen- 
values e n are given by 

' j (n + -) n = 0,l,... (3.12) 



The ground state energy E which is a sum over the first N energy levels 
clearly reproduces (|3.7f) . 

For later use let us also evaluate the density operator 

1 N - - 

p(0=vE^-O ( 3 - 13 ) 

JV i=l 



in the state (|3.8|) . The result is |10] 

Po,o(0 = i E^O 1 |0n(O| 2 = 2«7T-3/4[(l - )^ (1 _ ^ + o(1/JV)] 

(3.14) 

So far we have ignored the two-body interaction T-C 2 in (|3.6| ) (the addi- 
tional subscript on the ground state wavefunction and energy denotes this 
fact). We will see shortly that in the only sensible scaling available in the 
theory this term will be of lower order in N. In other words, the correc- 
tion terms to ( |3.7|) coming from the two-body interaction will turn out to be 
1/N lower order than N 3 ^ 2 , leading to the result that the total ground state 
energy is E = N 3 l 2 g /(2y/W) + o(iV 1 / 2 ). 

Collective field theory 

Since the result E oc iV 3//2 is rather unexpected (naively one would expect 
the pure Yang- Mills result N 2 ), let us try to understand this from some simple 
scaling analysis. The easiest framework to do such an analysis is collective 
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field theory ||11|| . The hamiltonian Ti in ( |3.6| ) corresponds to the following 
collective field theory hamiltonian (the subscript 'c' stands for collective) 



(3.15) 



n cYM = N\f\ /2*) ;^(p(0/2)[n 2 (e) + ^p 2 (0] 



Wei = N(gX /2n) Jd^p(0[2n^/\ + 1/2) 2 ] 
W c2 = N(fX /2n) (1/2Ao) / d£ J d£ K reg (£ - f /A )p(Op(f ) 

Here p(£) is the eigenvalue density defined in ( |3.13[ ). II(£) is defined such that 
n(Op(0 corresponds to the momentum density: / d£ II(£)p(£) = (1/N) J2 a Pa- 
This implies the following commutation relation 

[p(£),n(f)] = -JLdg5{t-?) (3.16) 

The theory is defined with the constraint 

dip(g) = l (3.17) 

which is a consequence of the definition ( |3.13j ). 
Classical analysis 

Let us regard the collective field hamiltonian TC C as a function of the 
classical variables p(£),IT(£). The Poisson bracket between them is simply 



(|3.16|) without the i on the right hand side 

m,m')}pB = -^dg8{Z-?) (3.18) 
This leads to the following equations of motion 

dm) = -^-[f + ^p 2 + f (f/Ao + l/2) 2 + (3.19) 

-^/df p(?)JW(£-?)A))] 
jvA 
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For time-independent solutions these reduce to 

¥Y p2 + 2 ~W m ° + 1/2)2 + W 2 I <*?P(?) ?)Ao)] = (3.20) 

If we take the naive N — > oo limit, the potential and interaction terms drop 
out and we get a constant density 

p(0 = constant (3-21) 

However, since that the range of £ is noncompact (due to the phenomenon 
of decompactification discussed in the last section) such a constant density is 
unnormalizable, that is, it cannot satisfy ( |3.17|) . There are no gauge- invariant 
cutoffs available on the range of £ either which can save the situation. 

Clearly in order to obtain normalizable solutions for p(£) the large N limit 
must be taken in such a way that one or both of the terms 7i cl and 7i c2 in 
( |3.15| ) are of the same order as H. c ym- Now it is easy to see that only 7i c2 by 



itself, being translation invariant, cannot produce localization. Indeed, if one 
includes in ( |3.2(J| ) terms coming from H c ym an d 7~tc2 an d drops those coming 
from Tici, the solution is still = constant, which is untenable. Thus we 
must ensure that the simple harmonic potential term in (|3.20|) survives in 



the large N limit (irrespective of what happens to the interaction term). In 
other words, in Q3.15|) 7i C YM an d H. c i must be of the same order. 

The clue to how the above can be achieved is provided by the formula 
( |3.14|) . This suggests that the correct coordinate in terms of which a sensible 



iV scaling of the collective field theory may be available is not £, but y, 
defined by 

x — y/Ny, equivalently ■=- + - = ay, (3.22) 

Ao 2 

where 

a = V 1/2 iV 1/4 = aVN (3.23) 
The new density variable p(y) is given by 

Kv) = P(M/dy = p{i)N l ' 2 a- 1 (3.24) 

Note that the density expectation value (|3.14|) in terms of the new variable 
reduces to 

PUV) = PoAON^a' 1 = 2vr- 3 / 4 (l - y 2 V^Y /2 0(l - y 2 ^) (3.25) 
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which is an iV-free expression. We also define a new momentum variable 
fi(y) by 

U(y) = n(OiV 1/2 a _1 (3.26) 

such that the Poisson bracket is kept invariant. This also ensures that iV 
scales out of the combination II 2 (£) + p 2 (£)/12 in H. c ym-> Eqm (|3.15| ). In 
terms of these, the collective field hamiltonian becomes 

n cYM = (g/2n) J dy (p(y) /2)[fl 2 (y) + (7r 2 /3)p%)] 

H cl = N^ 2 (g/2*) J dyp(y) 2th/ 2 (3.27) 

n c2 = NV* (g/27r)a 2 J dy f dy> p(y)p(y>)K Teg {a(y - y')) 

There are several remarks to be made here. First of all, note that TC c ym 
and Hci scale as iV 3 / 2 . It is easy to show that if we ignore the H c2 piece, 
then the time-independent equations (recall ( |3.20| )) in the new variables have 
a solution p(y) which exactly coincides with fl3.25|) and the the classical 
energy of this configuration is identical to (|3.7|). The second point is, 7i c2 
is subleading compared to the first two terms in 7i c by 1/N, if we demand 
that a is N- independent as iV — » oo. One can therefore treat this term 
as a perturbation. The reason for demanding iV-independence of a is that 
that is the only way a scaled form of 7i c2 can be obtained. Indeed one 
can easily rule out iV-dependence of a on other physical grounds also. If, 
for example, a grew with N, the function K Teg (a(y — y')) would become 
infinitely discontinous (its poles become infinitely dense). If, on the other 
hand, a went as some inverse power iV _7 ,7 > 0, then by using the formula 

oo 

Kr Cg (ay) = £ C(2m + l)(m/) 2m (3.28) 

771=1 

one can see that although from the point of view of convergence the T-C c2 
perturbation series would be sensible, it would have the unphysical feature 
that different parts of the two-body interaction contribute in different orders 
of perturbation theory leading to an incorrect representation of the nature 
of the interaction. 
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Scaling of length 

Eqn. ( p.23| ) is equivalent to 

L = — iV~ 1/2 3.29 
9 

In view of the fact that a is iV-independent, we see that as N —>■ oo, the 
(bare) length L of the space circle goes to zero as L = N~ l l 2 L where L is 



held constant. The parameter a can be identified as ygL/2%. 

Scaling of time and energy 

Let us rewrite the classical equations of motion using the variables p(y), 
U(y). We get 

N^d t p(y) = -(g/2n)d y (p(yMy)) 

N 1 / 2 d t U(y) = -(g/2n)d y [4- + ^p 2 + 2ny 2 + £ J dy> p(y')K rcg (a(y - y'))\ 

(3.30) 

It is clear that scales out of the classical equations if we work in terms of 
a scaled time i defined by 

t = N 1/2 i (3.31) 

Now a scaling of time (in the Euclidean framework a scaling of f3) implies an 
inverse scaling of energy 

E = N- 1/2 E (3.32) 

In other words, exp(iHt) = exp(iHt),exp(—/3H) = exp(—/3H). 
In this new scaling, the collective field hamitonian becomes 



n c = ^-J dy P -^[Ii 2 + ^p 2 +W + ^ J dy'p( y ')K Icg (a(y-y'))] (3.33) 
The scaled ground state energy is, therefore, 

E = iV 2 -^[l + o(l/iV)] (3.34) 

In the following we will study how the 1/N corrections follow from a WKB 
analysis. 
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1/N corrections and WKB 

The commutation relation between p(y) and II (y) is given by (recall 

\fi(y)My')]^-^9 y 5(y-y') (3.35) 

which implies that tl(y) can be represented by the following differential op- 
erator representation 

fi « = -w d *m (336) 

To construct a WKB solution we consider wave-functions of the form 

= exp(iN p S[p}) (3.37) 

where we have kept the real number p unspecified and will determine it by 
demanding a scaled time-independent Schrodinger equation. The latter is 
given by 

H c V[p] = m[p\ (3.38) 

where 7i c is given by (|3.33|) . Using the above differential operator representa- 
tion of Tl(y), ( p.36| ), we get the following equation determining S[p] in terms 
of E: 



£jdy'p(y')K Ieg (a(y-y'))} 



(3.39) 

In order to get a scaled Hamilton- Jacobi equation in the leading N order we 
must put p = 2. This implies the following perturbation series for S 

V[p] = exp[*iV 2 ,S], S = S + ^S X + -Ls 2 + ... (3.40) 

iv iv z 

It is easy to write down from ( |3.39 ) equations determining S n in terms of the 
lower order terms. 

The above analysis also implies that the energy E for the ground state 
should receive corrections to (13.71) as follows 



En = N 



9 1 A 1 A 

+ — A .i + — ^A , 2 + 



2^ N u ' N 2 



(3.41) 
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We have explicitly verified this by computing the many-body perturbation 
theory diagrams of the hamiltonian ( 3.15|) . 



Introduction of multiple flavours 

In the above we have considered the case of a single flavour, uf — 1- It 
is easy to extend the above calculations to multiple flavours. The essential 
new physics point is the following. The dirac sea \F) for np > 1 flavours of 
free fermions is given by 

ibfJF) = for n > n a J 

(3.42) 

iptJnlF) = for n < n a / 

where ?/>'s above represent the right-handed fermions. The left-handed fermions 
represent similar equations corresponding to the filling of all the levels start- 
ing from n a p + 1 upwards. 

If the fermi levels for different flavours are different, then one would an- 
ticipate several conceptual difficulties. For instance, if there is no notion 
of a single n F how does one construct £ a = A + n F X to describe the de- 
compactified effective theory of the eigenvalues? Fortunately, because of the 
chiral U(nf) x Uinf) symmetry of the hamiltonian, we must demand that 
the ground state is a flavour singlet. This happens only if the n a p l s above 
are all independent of i (this can be proved by operating the flavour rotation 
generators on \F)). Having said that, we can again use the notation \np) 
instead of \F) to represent the dirac sea. 

Using the above, it is easy to show that the dirac sea is again an eigenstate 
of the hamiltonian upto 1/N corrections. In other words ( |2.19|) is again valid, 
except that the energy eigenvalue Ep is now rif times the expression (|2.25|) . 
As a result in the calculation of the ground state energy according to (|3.6|) , 
Hcym remains the same but Ti. cl and 7i c2 get multiplied by rif. It is simple 
to show that the ground state energy is now given by 

E = N 2 [^W]-^= + o(l/ N)} (3.43) 



20F 

The eigenvalue density is given by 

p(y) = 2n 1 /\^\l - y 2 ^Wjfl 2 d{l - y 2 ^WJ) + (1) (3.44) 
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Eqn. ( |3.43| ) clearly shows that adding fermions affects the leading N be- 
haviour of the vacuum energy. 



3.2 Finite (3 partition function and Excitations 

The method presented above can be generalized to (3 ^ 0. In this case, 
one needs to consider excited states in the fermionic theory. This leads to 
subleading corrections (in 1/N) to the effective hamiltonian for the gauge 
fields. The free energy can again be computed. We will not present here 
the explicit result but merely state that the iV-scalings presented above also 
work at finite values of (3 and that the leading term (in N) in the expression 
for the free energy still has a non-trivial dependence on rif. 



3.3 Vacuum expectation value of the Wilson loop op- 
erator 

In order to understand how the introduction of dynamical fermions affects 
the leading large N behaviour of various quantities, it is instructive to calcu- 
late vacuum expectation value (W m ) = (TrU m ) of the Wilson loop operator 
(recall ( |1.8| ) and ( |1.9| )). Using the expression for the ground state wavefunc- 
tion Q2.29Q and the definition (|3.13 ) of the density operator it is easy to show 



that 

(W m ) = f dt(p(0)e i2 ^^° (3.45) 



On the right hand side the expectation value is meant in the ground state of 
the effective theory of eigenvalues defined by ( |3.5| ) and ( |3.6| ), and is given by 
(|3.44[) . Using (|3.22|) to express the exponent in terms of y we get 



r%2 

{W m ) = 2{-l) m (l + 7 —)J (x m )+o(l/N), Xm = 2n 3 / 4 a^ (3.46) 

where Jo denotes the Bessel function of order zero. 

Let us compare this with the result in pure Yang-Mills theory on a two- 
dimensional cylinder. The latter theory is described by the path integral 
( |1.3|) with no fermion integration and Sf = 0. It is easy to show in this case 
the eigenvalue density is a constant (which, unlike with fermions, is perfectly 
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allowed because the range of eigenvalues remains compact) and therefore 

(W m ) YM = S nfi + o(— ) (3.47) 

One can also recover this by taking rif — > in ( |3.46|) . 

It is clear that the introduction of fermions changes the result for {TrU m ) 
in the leading N order. 

4 Concluding Remarks 

To conclude we remark that QCD2 on the cylinder shows considerable amount 
of subtlety and surprise. The main surprise is that the leading term in a 1/N 
expansion for the vacuum energy and for the vev of the Wilson loop operator 
is changed by the existence of fermions, contrary to what one would expect 
by standard iV-counting arguments. The essential reason for this is that 
coupling to quarks leads to a 'decompactification' of the eigenvalues. This 
effect is not perturbative and persists even at N = 00, thus changing the 
leading behaviour of pure Yang-Mills theory. Indeed it is a very interesting 
question what the present analysis can tell us about 3 + 1 dimensional gauge 
theories in a compact space. In particular it would be interesting to know 
whether the phenonmenon of decompactification persists in four dimensions 
and whether as a result the pure Yang-Mills theory is affected by inclusion 
of quarks in the leading iV order. 
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